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ABSTRACT 

Time-resolved spectral analysis, though a very promising method to under- 
stand the emission mechanism of gamma-ray bursts (GRBs), is difficult to im- 
plement in practice because of poor statistics. We present a new method for 
pulse-wise time-resolved spectral study of the individual pulses of GRBs, using 
the fact that many spectral parameters are either constants or smooth functions 
of time. We use this method for the two pulses of GRB 081221, the brightest 
GRB with separable pulses. We choose, from the literature, a set of possible 
models which includes the Band model, blackbody with a power-law (BBPL), a 
collection of black bodies with a smoothly varying temperature profile, along with 
a power- law (mBBPL), and two blackbodies with a power- law (2BBPL). First, 
we perform time-resolved study to confirm the spectral parameter variations, and 
then construct the new model to perform a joint spectral fit. We find that any 
photospheric emission in terms of black bodies is required mainly in the rising 
parts of the pulses and the falling part can be adequately explained in terms of 
the Band model, with the low energy photon index within the regime of syn- 
chrotron model. Interestingly, we find that 2BBPL is comparable or sometimes 
even better, though marginally, than the Band model, in all episodes. Consistent 
results are also obtained for the brightest GRB of Fermi era — GRB 090618. We 
point out that the method is generic enough to test any spectral model with well 
defined parameter variations. 

Subject headings: gamma-ray burst: general — methods: data analysis — methods: 
observational 
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1. INTRODUCTION 

The spectrum of the prompt emission of a gamma-ray burst (GRB) is generally fitted 
with the celebrated Band spectral model (Band et al. 1993). This model represents a 
non-thermal spectrum, and it can be described in terms of two smoothly joined power 
laws. Generally speaking, Band model adequately fits most of the time-integrated prompt 
emission spectra of GRBs, though additional spectral components show up for some GRBs 
(e.g., Preece et al. 1996; Gonzalez et al. 2003; Shirasaki et al. 2008). But, to emphasise, 
the exceptions are very few in number in comparison with the large set of GRBs, which can 
be fitted with a simple Band only function (e.g., Kaneko et al. 2006; Nava et al. 2011). 
Zhang et al. (2011) have found that 15 out of 17 LAT detected GRBs could be fitted with 
Band only model (category I, in their notation). Hence, the Band function, till date, is the 
simplest, standalone model for GRB spectra, whether time-integrated or time-resolved. 

Although the Band model is statistically the most appropriate model of GRB data, 
its physical origin is yet to be identified. Over the years, many authors have investigated 
the underlying mechanism of the prompt emission. In the fireball model of GRB, the 
observed radiation during the prompt emission is attributed to a highly relativistic optically 
thick outflow, which thermalizes photons due to random collisions. This thermal energy is 
expected to be seen from the photosphere (Goodman 1986; Paczynski 1986; Meszaros & 
Rees 2000; Pe'er 2008), where the fireball becomes optically thin and radiation decouples 
from the matter. In this scenario, the light curve (LC) is expected to be a single, smooth 
pulse and the spectrum should be a Lorentz boosted BB, while the temperature of the 
BB adiabatically cools. Hence, the time-integrated spectrum should be a superposition 
of many BBs. Though some GRBs do have a single pulse, most GRB light curves are 
either superposition of many pulses or they are highly variable. Moreover, A BB in the 
Rayleigh- Jeans region has a photon index 1.0, which is much harder than that actually 
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observed for GRBs {a ~ —1). Many models have been proposed to overcome these 
difficulties. For example, the internal shock model (Rees & Meszaros 1994; Woods & 
Loeb 1995; Sari & Piran 1997; Kobayashi et al. 1997) assumes that the major radiation 
is not due to photospheric emission, but optically thin synchrotron radiation (SR) from 
internal shocks. One major problem with the synchrotron radiation is that the low energy 
photon index is limited to a < — | (Preece et al. 1998). Crider, Liang & Preece (1998), 
using time-resolved spectra of 99 GRBs, have shown that the instantaneous spectra and 
their evolution cannot be explained by SR — a often crosses the line of death, set by the 
synchrotron model. Another possibility is that the radiation is due to inverse Compton 
(IC) of the thermal photons near the photosphere (Thompson 1994; Pe'er et al. 2005, 2006; 
Beloborodov 2010; Lazzati & Begelman 2009). This process can indeed produce a Band 
like spectrum, but with a rather hard value of o; ~ 0.4 at the best, assuming slow heating. 

A unique prompt emission model of a GRB is yet to be settled. Prom the 
phenomenological point of view, the correct model can be identified by segregating the 
details from the average properties. For example, GRBs are superpositions of pulses 
(Nemiroff 2000; Norris 2005). Hence, one should use the individual pulses for spectral 
study, instead of the full GRB. The next step is to study the spectral evolution within 
the individual pulses. Hence, one should do time-resolved spectroscopy in order to extract 
greater information than merely an average spectral property of a pulse, e.g., average peak 
energy, isotropic energy etc. But, performing such a detailed study is difficult, as one loses 
photon counts. For example, Ghirlanda et al. (2010) have done time-resolved spectral 
study of 9 selected GRBs, detected by Gamma-ray Burst Monitor (GBM) onboard the 
Fermi satellite. As the photon count is low in each time bin, they could model the spectra 
only with a cut-off powerlaw. Parameters of a more complicated model, such as Band, 
cannot be well constrained from the time-resolved data. The solution to this problem 
lies in the realisation that spectral evolution is not totally unpredictable, and one can 
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suitably parametrize this evolution in order to reduce the number of free parameters of the 
description. For example, the spectral evolution of a GRB pulse is generally described as 
a hard-to-soft evolution (e.g., Liang & Kargatis 1996; Kocevski & Liang 2003; Nemiroff 
2012). Recently, Basak and Rao (2012a; b) have assumed this hard-to-soft evolution of the 
individual pulses of the set of 9 GRBs of Ghirlanda et al. (2010) to generate simultaneous 
spectral and timing model of the pulses, with essentially two parameters, namely, the peak 
energy at the start of the pulse {Epeakfi) and the characteristic evolution parameter (0o)- 
The basic assumptions in this approach, however, are not well established. For example, 
it was assumed that the spectral softening happens throughout the pulse, though there 
are evidences that some GRB pulses show a different behaviour like the intensity tracking 
spectral evolution (see e.g., Lu et al. 2012). Further, it was also assumed that the apphcable 
model is the Band function throughout. Hence, it is essential to critically examine all the 
applicable spectral models and their evolution to arrive at a correct pulse-wise description 
of a GRB. 

In this paper, we discuss a new method for pulse-wise spectral analysis where we 
parametrize the spectral evolution in order to arrive at the correct spectral description 
with a minimum set of free parameters. We apply this method to study GRB 081221, the 
brightest GRB with clean, separable pulses. We compare the results obtained for this GRB 
with those for GRB 090618 - the brightest GRB in the Fermi era. The organisation of 
this paper is as follows — in Section 2, we describe the analysis techniques and the basic 
assumptions of our model. Results are given in Section 3, and major conclusions are drawn 
in Section 4. 
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2. ANALYSIS METHOD 

2.1. Data selection and analysis 

The basic necessity for a good spectral analysis of GRB pulses is wide band coverage to 
identify additional spectral components. The Gamma-ray Burst Monitor (GEM) on board 
the Fermi satellite, with its wide band width and excellent sensitivity, provides a good data 
base for such studies. It has two scintillation detectors: the sodium iodide (Nal) detector is 
sensitive in the > 8 keV to ~ 900 keV range while the BGO energy range is ~ 200 keV to 
~ 40 MeV (Meegan et al. 2009). 

We examined the Nava catalog (Nava et al. 2011) of Fermi/ GBM GRBs and found 
that there are 112 bright (fiuence > 10~^ erg), long {5t > 15 s) GRBs and 11 of these GRBs 
have single/ separable pulses. GRB 081221 is the brightest among them. In Figure [H we 
have plotted the light curve (LC) of this GRB with Norris model (Norris et al. 2005) fitted 
for the two pulses. We have also made a systematic analysis of the other 10 GRBs and the 
time-integrated spectral analysis for all of them is given later. 

We use the CSPEC data for time-integrated study and the time tag event (TTE) 
data for the time-resolved spectral analysis. We choose 2 or more Nal detectors having 
high count rate and one/ both BGO detector (s). For source selection and background 
subtraction, we use the rmfit v3.3pr7 tool, developed by User Contributions of Fermi 
Science Support Center (FSSC). The background exposure time is chosen before and after 
the burst. This background is modelled by a polynomial of different degrees, according 
to the need. The PHA files are binned in energy channels so as to get a minimum count 
in each spectral bin. Typically, the Nal detectors are binned by minimum count > 40, 
while the BGO detectors are binned by requiring a minimum count of ~ 50 — 60. First, we 
perform time-integrated analysis for all the 11 bright, long GRBs having single/ separable 



-7- 



pulses. We implement both C-stat and minimization methods in rmfit. We fit either 
Band or Cut-ofF powerlaw (CPL) model. 

In Table 1, we report the best fit parameter values along with the corresponding 3a 
errors. The reduced C-stat and Xred '^i^h degrees of freedom (dof) are also reported. For 
comparison, we quote the results of Nava et al. (2011) for these GRBs. It is clearly seen 
that all these values are matching quite well. The sources of very minor deviations between 
the values of this work, done by C-stat minimisation and Nava et al. (2011) are: (i) 
mismatch between actual start and stop time, (ii) exact background selection and modelling 
and, (iii) the exact number of detectors and the channels used. Comparing the deviation 
of the parameter values, it is clear that the deviation due to using different statistics other 
than C-stat, i.e. minimisation, is much less than that due to these other reasons. Hence, 
we conclude that the statistics plays a minimal role in actual parameter estimation. In 
fact, the GRBs taken in our analysis are all bright GRBs (fluence > 10~^ erg). Hence, by 
default, the minimisation is a correct technique for parameter estimation of GRBs with 
high count rates. 



2.2. Spectral models for time-resolved study 

We select four models for the spectral study: a) Band model, b) blackbody with 
a power-law, c) a modified blackbody with a power-law and d) two blackbodies with a 
power-law. The Band function can be written in terms of the spectral indices {a and /3) 

and the peak energy (Epeak) as: 
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Here, Ab is the normalisation constant. A model consisting of thermal (blackbody with 
temperature kT) and non-thermal (power law with index F) components has been used 
earlier (see, e.g., Ryde 2004, Ryde et al. 2006, Ryde and Pe'er 2009, Pe'er and Ryde 2011). 
We name this function as BBPL. This can be written as: 

Ki X 8.0525£;2 ^ 

where Kl and K2 are normahsation constants. There are suggestions (e.g., Ryde et al. 
2010) in the literature of a modified blackbody (mBB), which may exist due to angular 
dependence of the optical depth and the observed temperature (Pe'er 2008). Hence, we 
also investigate this model with a powerlaw (mBBPL). The mBB model is a multi-colour 
blackbody disk model; the local disk temperature kT(r) is proportional to r~^. In several 
GRBs there are distinct additional thermal components (see eg. Shirasaki et al. 2008) and 
further, if the GRB spectrum is due to thermal inverse Compton (IC) of seed photons, 
then there may be, in principle, multiple photon baths. For brevity, we take one more BB 
component and call it 2BBPL model. This has essentially two BBs with two temperatures 
— kTfi and kTi and two normalisations — Kih and Ku. We use all the four models for our 
subsequent analysis. We note here that it is also possible to have combinations of the above 
models (like Band+PL, Band + BB + PL), but, at this juncture, we consider only these 
four generic models. This is essentially because of the fact that the time resolved spectra of 
GRBs generally consists of a broad peak with wings, which can be adequately captured by 
any of the above four models. 



2.3. Assumptions of the new pulse-wise spectral study 

Time-resolved spectral studies require a large number of parameters. If 'n' is the 
number of time bins, then a four-parameter model, such as Band or BBPL, requires 4n 
parameters for a full description. Our motivation, in this study, is to reduce the number 
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of parameters with some reasonable assumptions. The basic assumption we make is that 
the temperature (kT) and the peak energy (-Epeafc) follow smooth time evolution. This 
evolution has a break at the peak of the pulse (see also Ryde and Pe'er 2009). In the 
following discussion we shall refer to them as two episodes — rising and falling part. We 
assume that the time evolution law of kT and Epeak is simple powerlaw of time i.e., ~ t^, 
with different /i in different episodes. If 'm' and 'n' are the number of time bins in these 
two episodes, respectively, then this parametrization reduces the number of free parameters 
by 'm+n-2'. Note that, ideally one should use one more parameter to account for the start 
time, but we have chosen it to be zero for all the cases, except for the rising part of GRB 
090618, where the pulse starts from -1 sec. Hence, only for this case, we have chosen this 
start time to be -10 sec. As the parametrization is not corrected for the actual start time, 
the actual values of /i are not unique and hence, these values should not be used to compare 
different pulses. We also assume that for the Band model the photon indices (a and /3) are 
constants in each episode. Their values can be determined by simultaneously fitting all the 
time-resolved spectra in an episode, with the photon indices tied. For BBPL model, we 
assume that the index (F) of the powerlaw (PL) is constant in each episode. This index can 
be determined in the similar way as the photon indices of the Band. The parametrization 
of the norm of BBPL model is more complicated, as, unlike Band, we have two norms. In 
the Band model, norm is a free parameter. In order to have equal number of parameters as 
Band, either we can use an overall free norm with suitable parametrization of BBPL norms, 
or, we can parametrize one of the norms and treat the other as a free parameter. Ryde 
and Pe'er (2009) have shown that the parameter TZ = (Fbb / <^T'^Y^'^ either increases with 
time or remains constant. If we assume that the observed flux of the BB varies as a simple 
function of time as ~ t^, then the variation of TZ can be parametrized as ~ 
For fi < (/4, TZ will show the expected behaviors. The BB flux variation can be translated 
to the norm variation of Equation |2] as Ki ~ f^^ . Now, if we parametrize the BB norm 
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and treat the powerlaw (PL) norm as a free parameter of the model, we have the same 
number of free parameters as Band model. But, this makes the BBPL more constrained, in 
the sense that the overall norm is not a free parameter, as in the case of Band model. To 
overcome this difficulty, let us assume (which will be justified later) that the powerlaw (PL) 
component has norm variation as smooth as that of BB, i.e., K2 ~ f"^. Then the overall 
norm (K) can be made a free parameter by parametrizing the ratio of the BB norm and PL 
norm as K1/K2 ~ t". 

Thus, both Band and BBPL models have equal number of free parameters — 'm+n+8'. 
For Band model, 'm+n' are norms of Band function in 'm+n' time bins. The other 8 
parameters are a, (3, /i and the Epg^k at the starting bin, £'peafe(^o)) in the two episodes. 
Peak energy at any time 't' is determined by -Epeafc(^o) x {t/t^Y- For the BBPL model, 
'm+n' are overall norms (K) and the other four parameters are powerlaw index (F), /x, u 
and kT at the initial bin, kT(io), in the two episodes. For mBBPL model, we have an added 
parameter, namely, p, in each episode. Hence, the number of parameters is 'm+n+10'. For 
2BBPL model, the photons are boosted by the same material, hence, the BB parameters 
cannot be arbitrary. We assume that the ratio of temperatures and norms of these BBs 
are fixed, which should be determined by tying the ratios in all bins. Hence, the number of 
parameters in this model is 'm+n+12'. As an example, if there are 5 time bins preceding to 
the peak of a pulse and 10 bins afterwards, then the parametrization for Band and BBPL 
reduces the number of free parameters from 60 to 23; for a total of 25 bins this factor is 100 
to 33 and so on. 

To compare the significance of a model with respect to another, we have performed the 
F-test. The F value is defined in most general case as, F = ^l^f?^} , where index 1 is used 
for the primary model, while 2 is used for the alternative model. We compute the probably 
(p) of a given F value and thereby find the a significance (and the % confidence level) of the 
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alternative model preferred over the primary model. If the primary model is a subset of the 
alternative model, then the F value is defined as, F = an^j the difference 

between the dofs (M) is deemed as the dof of the primary model. 



3. RESULTS 

3.1. Time-resolved spectra of GRB 081221 

We start with a time-resolved spectral analysis so that some of the assumptions 
sketched above can be examined and validated. The major challenge in time-resolved 
spectroscopy is to define the time bin size, which crucially depends on two factors — a) 
timescale of spectral evolution and b) the minimum bin size allowed by the data, which in 
turn depends on the subjective decision of signal-to-noise ratio (SNR). One cannot violate 
the latter condition for a given SNR. The peak count rate of this GRB is ~ 4000, while we 
have demanded that each PHA bin should have at least 40 counts (SNR ~ 6.3, i.e., total 
~ 5000 counts, for 128 channels). Hence, we cannot choose smaller than ~ 1 sec time bin. 
First, we choose uniform time bins of 3 second to extract time-resolved spectra. Later, we 
reduce the time bin to 1 second to check any improvement due to finer time bins. 



3.1.1. Case I: Bin size of 3.0 seconds 

In Table HJ we report the results of Band and BBPL fit to the time-resolved data 
of 3 second bin size. The time bin starts from -1 sec, with 14 bins (numbered to 13). 
Approximately, first 4 bins belong to pulse 1, last 8 bins belong to pulse 2 and the 2 
intermediate bins belong to the overlapping region. For the BBPL fit, we first fit the 
spectra with the powerlaw index (F) free. We note that F is more or less constant for 
the major portion of the burst (bin #0-2 for pulse 1 and # 6-11 for pulse 2). We take 
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the average of F over these bins, separately for the two pulses, and found in both cases, 
r=1.83 with standard deviation (a) 0.14 and 0.10 for pulse 1 and 2, respectively. F has 
large error bars in the last bin of pulse 1 and last 2 bins of the second pulse. The values in 
the overlapping region (#4 and #5), apart from having large errors, may be ambiguous, 
and hence neglected. We freeze F to 1.83 for all the bins and redo the analysis. The 
corresponding values are also reported in Table HJ Note that, by doing this we are gaining 
one degree of freedom in each time bin. It is clear from Table H] that for the second pulse, 
the time-resolved spectral fit with the Band model are much better than those of the BBPL 
fit in terms of the Xred- ^'^^ contrary, these values are comparable in the first pulse. 
Hence, the first pulse may be dominated by thermal emission. Note, however, that x^g^ of 
BBPL is poor in bin of the first pulse, which has, in fact, the highest fiux. Hence, this 
points to a different radiation mechanism than a simple BB, which may appear in the high 
flux regions. 

In Figure [2l we have plotted the Xred BBPL with F thawed, BBPL with F frozen 
and Band flt with fllled circles, open circles and stars, respectively. The BBPL model is 
clearly inferior to the Band model for a major portion of the second pulse, specially in the 
regions where photon counts are high. In the flrst pulse (-1 to 11 sec), the BBPL model, 
except for the one bin is comparable to the Band model. We also flt mBBPL and 2BBPL, 
which are shown by fllled boxes and pluses. The mBBPL and 2BBPL are as good as the 
Band model, in terms of Xred- Hence, the correct model for the major portion of the burst 
is either of these three. 

3.1.2. Case II: bin of size 1.0 second 

In the time- resolved analysis of Section 3.1.1, the Band and other models show superior 
flts compared to the BBPL model in terms of Xred- But, as Zhang (2011) pointed out 
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for GRB 090902B, this can be an effect of evolution of the BBPL within a time bin. To 
investigate this, we make finer time bins of 1 second and redo the time-resolved spectroscopy. 
Note that this is the finest bin possible for a SNR ~ 6.3. In Table El we report the average 
values of Xred different model fittings, obtained for different bin sizes. Apparently, there 
are improvements in the Xred f*^^ finer bins, but we also see that these improvements are of 
the same orders for different models. Hence, it seems that the spectrum is not due to the 
evolution of BBPL, but one of these other models. But, all of these models are comparable. 
Hence, we cannot hope to find the correct model by taking finer bins. Presence of a simple 
blackbody (BB) is apparent in the first pulse, with the exception of the second bin, where 
one of the other models is correct. The second pulse is dominated either by a modified 
blackbody or a fully non-thermal radiation (Band) or inverse Compton (2BBPL). In order 
to find the right answer, one should carefully parametrize the spectral evolution in each 
episode of the pulses separately and reduce the set of free parameters in the description. 

3.2. Parameter evolution 

The spectral evolution in the pulses are not arbitrary. For example, the temperature of 
the blackbody evolves with time in a smooth way. This can be seen from Figure [21 where we 
have plotted the kT of BBPL (both F free and frozen) by circles. The peak energy variation 
is shown by pluses. The parameters evolve smoothly as a function of time. For a single 
pulse, Ryde and Pe'er (2009) have shown that the temperature remains constant or slowly 
declines with an average powerlaw index, (ax) = —0.07 with cr(aT) = 0.19, during the rise 
of the pulse and decays faster with an average index (bx) = —0.68 with cx^bT) = 0.24. The 
break time of this evolution has a strong positive correlation with the pulse peak time. We 
find a similar behaviour for both kT and Ep^ak evolution. Hence, the spectral evolution 
can be described by a simple time evolution of temperature {kT ~ t^) or the peak energy 
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{Epeak ~ t^)- The index, /i, in principle may have two values in the two episodes, namely, 
the rising and the falling part. 

Liang and Kargatis (1996; hereafter LK96) showed for FRED pulses that the peak 
energy of the EF(E) (or vF^j) spectrum follows a more complicated evolution. Ep^ak 
decreases exponentially with the running fluence as: 

TP /,\ 77 / <pBandit)\ , , 

Epeak[t) = Epeakfl exp (3) 

\ (pBandfl J 

where, Ep^akfl and (pBandfi are the constants of the Epeak evolution law, 
4'Band{t)= J^^ fit')dt', is the fluence at time t, f{t') being the flux. In Figure H] (up- 
per panels), we have plotted the In (Ep^ak) with the fluence. Note that the fluence of each 
pulse is calculated from their respective start time. Clearly, the Epeak evolution strictly 
follows the LK96 law in the first pulse. In the second pulse, however, the variation is not 
smooth throughout. In the falling part, the variation is clearly LK96 type, but in the 
rising part, the variation is rather "soft-to-hard". This effect may come due to overlap 
between two pulses. In fact, the first two bins of the second pulse belong to the overlapping 
region. Hence, they might be contaminated with the preceding pulse. However, the third 
bin, where there should not be any effect of the first pulse, also deviates from the LK96 
law. This might indicate that the second pulse is genuinely intensity tracking. Kocevski 
and Liang (2003) argued that the "intensity tracking" pulses for which this evolution does 
not appear very prominent, are rather made of more than one short hard-to-soft pulses. 
Ghirlanda et al. (2011), on the other hand, have analysed time-resolved spectra of 11 long 
and 12 short Fermi GRBs and found that the long GRBs appear to follow a "soft-hard-soft" 
trend, tracking the flux of the ORB, rather than a strict "hard-to-soft" evolution. Lu et al. 
(2012) have categorised GRB 081221 as one having a strict "hard-to-soft" pulse followed 
by "intensity-tracking" pulse. They have simulated overlapping pulses to show that in the 
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overlapping region, the spectral evolution may appear "intensity-tracking". However, they 
also found some single pulses to have "intensity-tracking" spectral evolution. Hence, the 
second pulse may be genuinely "intensity tracking". 

In Figure H] bottom panels, the temperature evolution is plotted against blackbody 
fiuence. Here the same behaviour is noticed. Hence, kT evolution of the first pulse and the 
falling part of the second pulse can as well be described by a similar exponential decay: 



t, f{t') being the flux at f. (pBBfl and kT^ are the constants of the evolution law. LK96 
law is empirical and a simpler version, in principle, can be used instead, such as a simple 
powerlaw of time. As fiuence is a monotonically increasing function of time, either of them 
can be used for evolution study. 

The BBPL model has two components. Hence, in order to parametrize the norms of 
this model one has to see the flux evolution of the individual components. In Figure [5l 
we have plotted both photon and energy flux of the individual components calculated for 
8-900 keV energy range. The flux evolutions look similar for both F free and frozen cases. 
Interestingly, the PL flux is as smooth as the BB flux, in each pulse. Hence, as argued in 
Section 2.3, we can safely assume that the ratio of their evolutions is a smooth function 
of time. In Figure El we have shown the evolution of a, (5 of Band and F of BBPL. The 
parameter /3, in many cases, have either large error bars or only an upper limit could 
be derived. In some cases, they peg to the value -10. It is clear from Figure [6] that the 
parameters remain reasonably constant at all episodes (rising and falling part) of a pulse. 
Hence, we tie them over all the time bins in a given episode to determine them with 
greater accuracy. This reduces the number of free parameters of the description of spectral 




(4) 



where (j)BB{t)=J^ f{t')dt' is the running fiuence of the blackbody component at time 
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evolution to a great extent, as described in Section 2.3. 

3.3. Results of the parametrized spectral fitting of GRB 081221 

The fact that the parameters are well behaved functions of time, makes the time- 
resolved spectroscopy more tractable, as we can reduce the number of free parameters in the 
pulse- wise description (see Section 2.3). Following the parametrization and tying scheme of 
Setion 2.3, we do the spectral analysis for the individual pulses of the GRB. In the following 
analysis, we use the TTE data of NaT — nO, nl, n2 and BGO — bO for our analysis. The 
constant, which takes care of the relative normalisation of the detectors should not vary 
throughout the burst. Hence, we freeze them to the values obtained in the time-integrated 
analysis, i.e., 2.25, 2.32, 2.34 and 3.24 respectively. Additionally, we make the following 
changes compared to the time resolved spectral analysis discussed earlier. We divide the 
data into spectra of equal total counts rather than equal time bins so that equal importance 
is given to all individual spectra. Further, the spectral data in each bin is regrouped into 
spectral channels to provide an uniform SNR. We also note that the 30 to 40 keV region 
of the spectrum of this GRB has the known calibration issues due to K-edge of Nal (see 
e.g., Guiriec et al. 2011). This does not matter much for parameter estimations, but, if one 
wants to compare different models in terms of ^ then it is wise to neglect these bins. In 
the following, we have done the spectroscopy by neglecting the 30-40 keV band. 

3.3.1. Analysis of Pulse 2 

This pulse constitutes the major portion of the burst. The count rate is > 3 times 
higher than pulse 1. Hence, we can analyse this pulse with greater accuracy and later use 
our experience to analyse the other one. We perform the analysis for two cases as follows: 
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Case I — Analysis for count per time bin > 3000: This analysis is done by dividing 
the second pulse from 17.0 s onwards, requiring > 3000 counts per time bin. We divide 
the pulse into two parts. 17 to 21.45 s is the rising part and the rest up to 40.45 s is the 
falling part. In the rising part, we get 3 time bins and the falling part has 9 of them. The 
spectral bins in the energies >100 keV sometimes show less than 2a count, while <15 keV 
show less than 3a count. Hence, we merge the 8 keV to 15 keV bins to form 1 bin; 100 to 
900 keV bins are merged into 7 bins, with progressively higher binning at higher energies . 
Similarly, spectral bins of the BGO (200 keV to 30 MeV) are merged into 5 large bins. All 
the spectral fit parameters are listed in Table HI 

(a) Rising part: In the rising part (first four rows of Table Hj), the BBPL, compared 
to the other models is inferior with (dof) — 455.92 (354). If we parametrize only the 
BE norm, and treat the PL norm as free parameter, the corresponding (dof) is 487.68 
(354). A BBPL fitting with no constraint gives (dof) = 451.09 (348). While we have 
gained dof, the Xred remains of the same order, which confirms that the parametrization 
works. In comparison to BBPL, mBBPL is a better fit with (dof) = 355.68 (353), and 
significance of 2.55 a (98.93% confidence). Band is better than BBPL model, with x"^ (dof) 
= 364.41 (354), and significance of 2.37 a (98.23% confidence). This suggests that the 
radiation mechanism in the rising part of the second pulse may be a photospheric emission, 
but the thermal part is a modified blackbody (mBB) rather than a simple BB. However, if 
we compare these values with a 2BBPL model, then we immediately see that this model is 
the best with x^ (dof) = 351.78 (352). As the set of parameters of BBPL model is a subset 
of 2BBPL, the significance of 2BBPL compared to BBPL is much higher — 9.29 a (100% 
confidence). Compared to the Band model, 2BBPL model has a significance of 0.86 a 
(60.94% confidence), which shows that 2BBPL is only marginally better than Band model. 

(b) Falling part: In the falling part (see Table H]), the Band model is better compared 
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to the BBPL as well as the mBBPL model. Compared to the mBBPL model, the Band 
model has 35.22 less with one more dof. A comparison with 2BBPL model, on the 
other hand, shows that 2BBPL is better than Band model at 1.03 cr significance (69.71% 
confidence). Compared to Band model, 2BBPL has 40.64 less with two less dof. The 
Band model does not show much difference in terms of residuals of individual spectral fit. 
But, it is only when we perform a parametrized joint fit, then we realise that the 2BBPL 
model is marginally better than Band model (Table 4). Hence, in this region either Band 
or 2BBPL is the best model, with 2BBPL marginally better than Band model. In Figure [TJ 
we have shown the significance of 2BBPL fitting over the BBPL fitting study. 
The residual of BBPL model shows excess at various channels. No such structure is visible 
in the residual of 2BBPL fit. Note that the Nal K-edge is present in both the residuals 
between 30-40 keV. We have done the fitting both by including and excluding this band. 
When all the channels are used the (dof) of BBPL and 2BBPL are 340.85 (217) and 
239.82 (215) respectively. 2BBPL is preferred over BBPL at a significance of 8.42 a (100% 
confidence, p=3.86 x 10^^^). If we exclude the 30-40 keV bins, the corresponding x^ (dof) 
are 300.26 (197) and 193.17 (195), while 2BBPL is preferred at a significance of 9.01 a 
(100% confidence, p=2.10 x lO^^^) 

Case II — Analysis for count per time bin > 1000: To check whether lowering the size 
of the time bins improves the BBPL fitting, we perform the same analysis for count per 
time bin > 1000. As the count rate is lower, we merge the 100 to 900 keV of Nal detectors 
into 5 channels rather than 7. The rest of the binning remains the same. We obtain 10 bins 
in the rising and 29 bins in the falling part. 

(a) Rising part: The values of x^ (dof) for BBPL, Band, mBBPL and 2BBPL are 
1247.91 (1187), 1328.50 (1187), 1239.20 (1186) and 1222.76 (1185), respectively. Compared 
to BBPL model, the Band model is preferred at 1.36 cr (82.76% confidence), mBBPL 
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is preferred at 1.56 cr (88.17% confidence), while 2BBPL is preferred at 9.66 o (100% 
confidence). 2BBPL is preferred over Band model at 0.89 o (62.60% confidence). Hence, 
the conclusions of Case I remains unaltered. 

(b) Falling part: Similarly, in the falling part, the finer bin makes equal impact on all 
the models and hence, the conclusion remains the same. Note that compared to mBBPL, 
2BBPL model has 113.84 less x^, with the cost of one more dof. Hence, 2BBPL model is 
preferred over mBBPL in the falling part (1.31 a with 68.99% confidence). Compared to 
BBPL model. Band, mBBPL and 2BBPL are preferred at 3.12 a (99.82% confidence), 2.67 
cr (99.24% confidence) and 19.61 cr (100% confidence), respectively. 2BBPL is marginally 
better than Band at 0.95 cr (65.64%). 

3.3.2. Analysis of Pulse 1 

The time-resolved spectra of this pulse are extracted by requiring > 1000 counts per 
bin, as the photon count is ~ l/3rd of pulse 2. As before the spectral bins of Nal are 
binned in 8-15 keV and 100-900 keV, while BGO spectral channels are merged to form 5 
broad channels. The 30-40 keV band is neglected. The results of different fits are reported 
in Table [51 

(a) Rising part: The rising part has only one bin from -1.0 to 2.15 s. Interestingly the 
BBPL model is marginally better than the Band model in the rising part (see Table at 
0.87 cr (61.46% confidence). The mBBPL model has comparable as BBPL, but with one 
more parameter. 2BBP1 model has a significance of 1.04 cr (70.17% confidence) compared 
to Band model, while the same model has a significance of 2.19 a (97.12% confidence) 
compared to BBPL model. 

(b) Falling part: In the falling part of this pulse, the BBPL model is no longer the 
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best model. The Band is the best model with (dof) = 544.65 (473). mBBPL and 
2BBPL models are comparable to Band with ^ (dof) 548.22 (472) and 544.79 (471). 
Hence, the spectrum may be still thermal, though the thermal part is no longer a simple 
BB, but either a multi-color BB (mBB) or has multiple spectral component (one more BB) 
or simply synchrotron dominated (Band). Note also that the low energy photon index 
[a = — 0.86^Q^g), in the falling part is within the regime of synchrotron model, which is 
clearly in contrast with the rising part (a = — 0.55lo22)- This phenomenon of softening of 
photon index at the falling part of a pulse can be seen for all the pulses (see Table H] and 
IS]). In Table ini we have listed all the significance levels (in terms of p- value, sigma level 
and confidence level — CL.) of a model over another. Modeli is the primary model, while 
Model2 is the alternative model. It is clear from the table that 2BBPL model is preferred 
over Band model, though marginally, in some cases. The p values denote the probably that 
the alternative hypothesis is incorrect. Hence, lower the value of p, better is the alternative 
model over the primary model. Only in one case, namely the falling part of pulse 1, Band is 
preferred over 2BBPL. But the p- value of this case is 0.48, which signifies that they are only 
comparable. Interestingly, if we use finer bin size, the significance of Band and mBBPL over 
BBPL model decreases in the second pulse. The significance of 2BBPL, however, increases. 



3.3.3. Connection between the rising and the falling part 

The smooth variations of the parameters demand that the temperature, kT (of BBPL 
or mBBPL or 2BBPL) or peak energy, Epeak (of Band model) should be continuous 
functions of time, even during the pulse peak time. Hence, these parameters should match 
at the peak within errors to the one predicted by the empirical law. We follow the evolution 
law of the rising part to predict these parameters in the first bin of the falling part. We 
compare these values with the corresponding observed values. In Figure [H we have plotted 
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the observed values with respect to the predicted values. Note that the error bars of the 
observed values are much less compared to those of the predicted values. The sources of 
errors in the predicted values are errors in the evolution parameter, /j, and the errors in 
the actual parameter at the starting bin of the rising part. Generally, the parameter /i has 
large errors, which affects the errors of the predicted values considerably. The data points 
are essentially the same for both wider bin (open circles) and finer bin (filled circles). The 
dot-dashed line, which shows the equality of the observation and prediction, goes through 
all the points. 



3.4. Comparison with GRB 090618 

GRB 090618 is an interesting GRB in many aspects (for details see Ghirlanda et al. 
2010; Rao et al. 2011; Basak & Rao 2012a). This is a GRB with very high fiuence. The 
redshift (z) of this object is 0.54 and the total fiuence is 3398.1 ±62.0 x 10"'' erg cm~^, when 
integrated over its duration (182.27 s). In terms of fiuence, this is the brightest among all 
GRBs detected by Fermi. This is a very long GRB with multiple peaks. It has four broad 
pulses as follows. Pulse 1 is rather a clean precursor from -1 s to 40 s. The second pulse is 
well separated from this precursor, and occurs from 50 s to 75 s, with two structures in 50 
to 61 s. The third pulse, occurring from 75 s to 100 s, is contaminated with the faUing part 
of the second pulse, and the rising part of the fourth pulse. The fourth pulse occurs from 
100 s to 124 s. Though the secondary pulses (i.e., other than the precursor) are sometimes 
overlapping, we can still examine the spectral variation in the first pulse and in the major 
portions of the other pulses. GRB 081221 has only one secondary pulse, which makes it 
more convenient. However, in contrast with GRB 090618, the precursor of GRB 081221 has 
overlaps with the secondary pulse. In order to compare the results of GRB 081221, we shall 
take the precursor and the second pulse of GRB 090618. 
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3. 4-1- Precursor pulse 

The time- resolved spectra of this pulse are obtained by requiring minimum of 1000 
counts per bin. We obtain 10 spectra in the rising part (-1.0 to 14.15 s) and 11 spectra in 
the falling part (14.15 to 40.85 s). For the rising part, the parameters, fi and u are obtained 
by assuming the start time at -10.0 s (see Section 2.3). The results of spectral fitting by 
different models are reported in Table [71 The advantage of this pulse 1 over 081221 is it 
is longer and brighter, enabling us to parametrize the rising part. Also, this pulse is fully 
separated from the secondary events. It is clear from Table [7] that BBPL model is inferior 
to the mBBPL model for this pulse. In case of GRB 081221, we found that the BBPL and 
mBBPL models are comparable to each other, and marginally better than Band model. In 
this case, we definitely need a mBBPL rather than a BBPL for comparable fit as Band. 
Note that 2BBPL model is only comparable, but not better than mBBPL. Hence, the 
spectrum in the rising part may be modified blackbody dominated. In the falling part, 
mBBPL is again comparable to Band model. The 2BBPL model is superior than all the 
models in the falling part. The same conclusion was drawn for GRB 081221. Hence, there is 
hardly any difference of spectral evolution in the precursor pulse between these two GRBs. 

3.4-2. The second pulse 

This pulse is more difficult to analyse, because, it has two small structures in the 
rising part and two pulses at the peak. Hence, we ignore the 50 to 61 s of data and analyse 
only 61 to 75 s. The falling part covers 64.35 to 74.95 s, where we obtain 24 time-resolved 
spectra requiring 2000 counts per bin. This, in principle, can be lowered as we have higher 
count rate, but, following the analysis of GRB 081221, we restrict ourselves to a moderate 
count per bin. The results of joint spectral fit with various models are reported in Table |Hi 
It is clear from this table that BBPL and mBBPL are not the correct models in the falling 
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part. In comparison Band model is much better. However, the 2BBPL model, which shows 
comparable or better fit than Band in all cases, is only comparable to mBBPL in this 
particular case. This may arise due to the fact that this pulse is actually a combination 
of two highly overlapping pulses (see Rao et al. 2011). Hence, we redo the analysis on 11 
spectra from 69.25 s to 74.95 s, which covers only the falling part of the second pulse. We 
obtain the following (dof): 1137.67 (991), 1179.13 (990) and 1156.25 (989) for Band, 
mBBPL and 2BBPL, respectively. Hence, the Band and 2BBPL models are comparable, 
though we cannot rule out the possibility of contamination even in this falling part. 

In the rising part (see Table [H]), though mBBPL is better than BBPL, it is inferior to 
2BBPL and Band. Hence, the rising part of this pulse is probably synchrotron dominated. 
2BBPL is comparable to Band. In summary of this pulse, the whole episode can be 
described by Band model. This is in contrast with the second pulse of GRB 081221, where 
mBBPL clearly dominates the rising part, and then it is taken over by Band. 



4. Discussion and Conclusions 

From purely data analysis point of view, there are essentially two parts in a GRB 
spectrum — the peak of the spectrum and the wings, which extend to very low and very 
high energies. The fundamental difference between a Band only model and the models 
with thermal and non-thermal parts is: while Band accounts for the peak position, with 
exponential fall in the wings, the other models have thermal component to account for the 
peak position, and a powerlaw, falling slower than an exponent, holding the spectra at the 
wings. In principle, all these spectra should be equally good at the peak position, except for 
the fact that Band and mBB have broader peak than simple BB, while 2BBPL has double 
hump. Hence, the difference between these models arises mainly in the wings. The BB is 
inferior to the others if the peak is not narrow. As photon count at the peak is larger, the 
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residual should show up immediately. But, this is not easy to see if the difference occurs 
in the wings. Consequently, the three very different models, namely. Band, mBBPL and 
2BBPL show comparable x^, while fitting time-resolved data. Hence, re-binning at these 
wings plays a very important role to pin down the correct model. However, we cannot 
expect an order of magnitude improvement in the x^, because, binning in these wings gives 
6-8 broad bins (see Section 3.3.1) with large errors, while the major statistics comes from 
the peak position. 

We have found, in our analysis, that the spectrum changes from one pulse to the 
other, and even within a pulse. The fact that one of these four models is superior than the 
others, in a particular episode, points to the fundamental radiation mechanism. We see that 
this change of superiority is not random. For example, the first pulse of both the GRBs 
have shown that a mBBPL model is better, though marginally, than Band in the rising 
part. For GRB 081221, this could be described even by BBPL model in the rising part, 
which is really pointing towards the thermal origin of radiation in the first pulses. Similar 
observations are reported in the literature e.g., Ryde et al. (2009) showed, in the 1-3 s time 
bin of BATSE detected GRB 981021, that a BBPL is better fit than a Band model. Note 
that the low energy photon index (a) of the Band model crosses the synchrotron limit in 
the rising part, where the thermal models are adequate. However, in the falling part of all 
the pulses, where Band is better than mBBPL, a is consistent with the synchrotron limit. 
Hence, we can safely conclude that the radiation mechanism starts with a thermal origin, 
but is rapidly overtaken by synchrotron mechanism. The first pulse may be dominated 
by the photospheric emission in all episode, but the second pulse is mostly synchrotron 
dominated. The second pulse may or may not have a thermal origin. For example, the 
second pulse of GRB 081221 shows a mBBPL model in the rising part, which then becomes 
synchrotron dominated in the falling part. On the other hand, the second pulse of 090618 
is always synchrotron dominated. Hence, the transition between these different radiation 
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paradigms is smooth and repeatable. 

In comparison to the mBBPL model, we note that the 2BBPL is particularly better in 
all episodes. This model sometimes shows superiority to the Band model, even at the falling 
part of a pulse, except for the second pulse of GRB 090618, though we cannot rule out 
the possibility of two highly overlapping pulses in this particular case. Softer component 
than Band were reported for a few BATSE GRBs by Preece et al. (1996). Shirasaki et 
al. (2008), using the time-resolved spectral data of GRB 041006, detected by HETE-2 (2 
keV to 400 keV), found multiple spectral components, each having characteristic evolution. 
After the launch of Fermi satellite, these earlier claims were reconfirmed in some cases. 
For example, Guiriec et al. (2011), fitting the time-integrated spectrum of GRB 100724B, 
have shown the presence of an additional blackbody (BB) component along with the 
traditional Band spectrum (also see Burgess et al. 2011). In our analysis, we have used 
two blackbodies to account for the softer components. The origin of these two components 
is speculative. They might be different locations of the boosted front of the fireball having 
same temperature, but different boosting factors. Alternatively, they can be different seed 
photon baths, up-scattered by the bulk material. Irrespective of its origin, this model shows 
superiority to all other models in all episodes. Note that, though 2BBPL model has double 
hump in the peak, one of the peaks may occur in the lower wing (i.e., <15 keV). Hence, 
it is easy to identify this model, only if the difference occurs at the peak. Figure [7] clearly 
shows the double hump in the residuals of the BBPL fit. Hence, it is easy to visualise the 
2BBPL model from this figure. The Band model, however, has similar residuals as the 
2BBPL model. Hence, the data are not sufficient to distinguish between these two models, 
except when we perform a parametrized joint fit. In Figure [9l we have shown the marginal 
superiority of the 2BBPL fit over the Band model study of the rising part of the 

first pulse (-1.0 to 2.15 s) of GRB 081221. In the right panels, we have plotted the fitted 
data neglecting the 30-40 keV channels. The upper panels are 2BBPL fits, while the lower 
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panels are Band model fits. Residuals of the Band model shows structures with excesses 
in 15 keV, 50-60 keV, and 150 keV regions. These are not present in the residuals of the 
2BBPL model. Of course, the difference is not as prominent as the case of Figure [71 The 
2BBPL model is preferred over the Band model at 1.04 a (p=0.298, 70.17% confidence). 
Hence, 2BBPL is only marginally better than the Band model. 

To visualise the evolution of the lower black body component, we have plotted in 
Figure [TOl the residuals of 2BBPL fit, with the lower BB omitted. This technique is well 
known for finding iron line profile in the inner accretion disks of black holes (see e.g.. 
Miller 2007). We fit the spectrum with the 2BBPL model and then omit the lower BB. 
The residual (expressed as normalized counts keV^^ s^^) of the fit clearly shows this BB 
component. Residuals of different detectors are shown by different markers. We have 
overplotted the lower BB model (in terms of normalized counts keV^^ s^^) on the residual 
to guide the eye. We have plotted these residuals for second, sixth and ninth time bins from 
top to bottom panels to show that the BB peaks at lower energies at later times. 

In summary, we have rigorously used the evolution of parameters in the pulses of a 
GRB to construct various spectral models with a minimum number of parameters. We 
have constructed Band model with parametrized peak evolution and tied photon indices, 
BBPL with parametrized norm ratio of the BB and PL, parametrized temperature, and 
tied PL indices. Apart from these, we have used mBBPL and 2BBPL, which, other than 
the parametrizations of BBPL, have tied p indices, and tied ratios of temperatures and 
norms, respectively. This new method is quite general in the sense that any such model 
can be incorporated with suitable parametrization. The fact that the parametrization 
works demands a close look into the theoretical predictions of various radiation models. 
These models, irrespective of their complexities, should produce such smooth variations of 
parameters within a pulse of a GRB. Also, if there is really a transition from one radiation 
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mechanism to another, one should correctly model the mechanism of such transition. The 
fact that the synchrotron model is applicable at the falling part of the pulses, without 
invoking any other component, is really intriguing and demands a close look at the 
predictions of the internal shock model. 

One of the surprising result obtained in this work is that the 2BBPL model is 
statistically superior to the other models in most of the episodes in these two GRBs. Basak 
& Rao (2012c) have used this model for GRB 090902B. The residual of BBPL fit clearly 
shows double humps (see Figure 3 of Basak & Rao 2012c) which are taken care by the 
two peaks of the 2BBPL model. We selected this model purely in a phenomenological and 
data analysis perspective: distinct blackbody components (apart from the main peak in 
the spectrum) are seen in a few GRBs and while looking at the residuals, two humps are 
clearly discernible in a few time bins. Since these two are the brightest GRBs for such 
analysis (GRB 081221 is the brightest GRB in the category of GRBs with single/ separable 
pulses and GRB 090618 is the brightest GRB in the Fermi era), it is unhkely that we can 
reinforce this result by analysing data from other GRBs. One method could be to get the 
pulse-wise spectral parameters of a sample of GRBs and relate them to other properties of 
GRBs like redshift, afterglow properties etc. This will not only help us to identify the most 
appropriate spectral description but also to identify the emission mechanism operating 
during the prompt emission. 
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Table 1:: Results of time-integrated spectral analysis of 
the GRBs. 



GRB 
(Model) 


^1, ^2 


This work 


Nava et al. (2011) 


G-stat 




080904 
(GPL) 


-4.096, 21.504 


a=-1.22l°:2i 
Ep = AQ.lt^^ 
&=l 08 f597l 


-^p = 39.8i3;34 
^=1 23 f597l 


a=-1.14±0.05" 
Ep = 39.24 ± 0.75 
G=l. 14(587) 


080925 
(Band) 


-3.840, 32.0 


(5 ==-2.34+?:5° 

Ep — 158.9^24^4 
G=l 17 f712) 


/3 ==-2.24^^71 
Ep — 157.3l24!9 
^=1 14 f712) 


q;=-1.03±0.03 
13 ==-2.29±0.08 
Ep = 156.8 ± 7.07 
G=l. 13(716) 


081118 
(Band) 


0.003, 19.968 


«=-o.42i°:l^ 

/3 ==-2.18l°i^ 
E„ = 55 93+?o? 

j-^p '-"-'•^'-'—12.5 

C=1.17 (716) 


a=-0.37;°:I° 
15 ==-2.141°;!^ 

E„ = 54 0"^!o n 

xLrf=1.02 (716) 


a=-0.46±0.10 
/3 ==-2.29±0.05 
Ep = 56.79 ± 2.77 

C=l. 16(601) 


081207 
(Band) 


0.003, 103.426 


a=-0.58+°-;,° 
/3 ==-2.151q;33 
Ep = 363.4^™;^ 

G=1.43 (713) 


a=-0.58l°;i? 
^ ==-2.13l°i? 
Ep = 364.5;ii 
xL=l-02 (713) 


a=-0.58±0.02 
p ==-2.22±0.7 
Ep = 375.1 ± 13.2 
G=l. 74(596) 


081217 
(GPL) 


-28.672, 29.696 


a=-1.09l°:i^ 
£'p = 193.0l;37;3 

C=1.19 (715) 


Q;=-1.10l°;i^ 
Ep = 200.5^11:? 

xL=l-06 (715) 


q;=-1.05±0.04 
Ep = 189.7 ± 11.2 

C=1.46(599) 


081221 
(Band) 


0.003, 39.425 


«=-0.84lS:^^ 
^ ==-4.24tV^ 


«=-0.84lH^ 
^ ==-3.89l?:f 


a=-0.82±0.01 
p ==-3.73±0.20 



Continued on next page 
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Table 1:: Results of time-integrated spectral analysis of 
the GRBs. 







Ep — 85.25l3;Q8 
G=1.64 (595) 


Ep = 85.0911^3 
.=1.49 (595) 


Ep = 85.86 ± 0.74 
G=1.67(600) 


081222 

(Band) 


-0.768, 20.736 


q;=-0.891°:}^ 

/5 ==-2.46l?:i? 
E„ = 169.2+^H 
G=1.12 (595) 


a=-0.89l°i^ 
/3 ==-2.32;°^^ 

En — 168.9^00 o 

X?,,=1.07 (595) 


a=-0.90±0.03 
P ==-2.33±0.10 
Ep = 167.2 ± 8.28 

G=1.23(604) 


090129 
(GPL) 


-0.256, 16.128 


a=-1.43l°:i^ 
E, = 170.4^ 
G=1.09 (596) 


a=-1.46t°:}^ 
Ep = 195.51^^2^ 
.=1.03 (596) 


a=-1.46±0.04 
Ep = 166.0 ± 15.1 
G=l. 12(602) 


090709 
(GPL) 


0.003, 18.432 


«=-i.04iS:i 

Ep = U6.7t'^i 
G=1.05 (596) 


a=-1.08l°J} 
= UA.ltlZ 
X,^ed=l-01 (596) 


a=-0.96±0.08 
Ep = 137.5 ± 12.5 
G=l. 17(602) 


091020 

(GPL)'= 


-3.584, 25.088 


a=-1.31+°-?^ 

— U.lo 

Ep = 255.7^92.0 
G=1.03 (354) 


1 OO+0.22 

— u.iy 
= 276.4lf«?:° 

xL=0-95 (354) 


a=-1.20±0.06 
/3 ==-2.29±0.18 
Ep = 186.8 ± 24.8 
G=l. 18(354) 


091221 
(Band) 


-2.048, 37.889 


/3 ==-2.40l°i° 

Ep — 191.3^47 5 

G=1.42 (474) 


«=-0.62l°i^ 

/3 ==-2.26tltl 
Ep = 189.51^?:^ 

XL=1-12 (474) 


a=-0.57±0.05 
P ==-2.22±0.10 
Ep = 194.9 ± 11.6 

G=1.44(466) 



" The errors quoted from Nava et al. (2011) are symmetric errors. Errors for this work are 3a errors. C is the reduced C-stat value, the 
number in the parentheses axe dof. The Band spectrum showed unbound 3a errors, we found better fit with CPL for this GRB 



Table 2: Results of BBPL and Band fitting of time-resolved data of GRB 081221. For BBPL, values are quoted for both 
powerlaw index (F) free and frozen to mean value at the high count rate region — 1.83. Ki is BE normalization, while 
K2 is that of the PL. The norms should be used as the relative normalization, as there is another constant multiplication 
due to detector effective area. Bin ^^s are 0-13, with denoting -1 to 2 s and subsequently equal bin size of 3 s is 
applied. Errors in K2 of the bin 7^ 13 of BBPL (F free) could not be determined. In many cases, only the upper error 
in /3 could be determined. 



Bin 



BBPL (r free) 



BBPL (F frozen to 1.83) 



Band 



kT 



kT 



Ki 



,(<lof) 





1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 



03 



+4.91 
4.28 



.26 



+2.45 



14 



82 



+3.12 
2.07 

a+2.56 
'-2. 14 

+1.69 
1.13 



36 



+2.01 
-1.74 



6V 



69 



+0.74 



77" 



41 



+0.84 
0.82 

+0.81 



+ 1.17 
■1.05 

-1.19 
-1.06 

-1.19 
-1.07 

+1.36 
1.11 



2.99 

1.76 
0.78 
0.74 
0.61 



+0.67 
-0.68 

+0.39 

0.36 

+0.28 
0.27 

+0.31 



+0.28 
0.26 



1 -19+0.30 
^■^^-0.29 



8.52 



+0.68 



13.67 



+0.77 



8.93 
5.93j 



+0.60 
0.58 



3.36 



+0.38 
0.37 



1 Q7+0.31 
-0.30 

. 42+0-29 
^•^^-0.27 



0.65 



+0.17 
0.19 



1.73 

1.77! 
1.99 
2.16 
1.85 
2.08 



+0.42 
-0.25 



+0.27 
0.23 

+0.40 



+0.27 
0.34 

+0.16 
0.13 



1 77+0.07 
' -0.06 



1.73 
1.761 
1.86j 
1.84 



-0.06 
-0.05 

-0.05 
-0.04 



-0.06 
-0.06 



2 01+°-l* 
^•L'l_0.12 

9 91+0.36 
^■^-■■-0.24 



3.65 



+3.10 
1.58 



5.31 



-12.72 
-3.15 



14.27^ 



20.66 
32.99^ 
9.05: 
49.09; 



28.87 
-13.36 



-17.99 
-7.42 



+31.10 
-19.27 



35.88l?:«^ 



41.65] 
53. 13^ 



6.89 
9.22 



75.42tli;^« 



54.40 



+ 13.64 
11.52 



40.27t?2.8i 
37 04+^^ °* 
164.92 



1.03(67) 
1.34(76) 
1.02(67) 
0.94(69) 
0.92(134) 
0.91(147) 
1.56(178) 
1.86(182) 
1.87(180) 
1.69(175) 
1.58(164) 
1.32(149) 
0.99(141) 
1.13(127) 



38.28 
16.92 

9.07 
9.07 
6.76 
9.18 



+4.78 
-3.97 

-1.87 
-1.63 

+1.49 

-1.33 

+ 1.39 
1.32 

+ 1.17 
1.08 

+1.06 
1.02 



23.10l°;«2 

24.341^5^ 

20.51±!lir 



14.06Z 



12.41 



+0.78 
-0.72 



10.27l0;^3 



8.96; 



+0.74 
0.72 



7 90+1.02 



, 14+0.49 
•3-J^*-0.46 



1.83 
0.84 
0.80 
0.63; 
1.19 



+0.34 
-0.33 

-0.26 
-0.25 

+0.26 
0.24 

+0.21 
0.20 

+0.27 



8.92 

14.53 
9.46 



-0.26 
+0.60 



0.49 

0.58 



+0.50 
0.49 



5.87Z 



3.35 



+0.37 

0.37 



9 r,Q + 0.30 
Z.UO_Q 29 



1.50: 



+0.26 
-0.25 



0.65 



+0.20 
0.20 



7.18 



+2.25 
2.16 



17.33^ 



10.94 
9.21 
8.13 



+3.28 
3.10 

+3.02 
2.87 



+2.89 
2.71 



18.14t^;J? 
43.4713-37 

+3.57 



55.80 



66.35 
67.70l 



+3.95 



51.03 



-3.92 
-3.86 



8.31 



+2.91 
2.79 



n 44+2-39 



1.01(68) 
1.33(77) 
1.02(68) 
0.99(70) 
0.91(135) 
0.97(148) 
1.56(179) 
1.89(183) 
1.90(181) 
1.69(176) 
1.57(165) 
1.36(150) 
1.03(142) 
1.15(128) 



-0.28t0;|« 



-0.69^ 



-0.24 

-0.85j 



+ 1.64 
1.10 



n KO+3.16 
^■^^-1.60 

-1 06+°-^^ 
l.UO_g 32 



+0.09 
-0.09 

-0.08 
-0.08 

-0.09 



-0.45_ 

-0.3lj 
-0.61^ 
-0.88^ 
-1.02 



+0.11 
0.10 



-0 61+°-*^ 
L'-"Jl-0.31 

"■^-'■-0.37 

1.18+M? 



-10.0 



-3.76: 
-2.55i 



-1.30 

-00 

-0.29 



_292+o^46 

-10.0 

-3.821^^^^ 

•^-■J^-1.09 

_937+i9.37 

_g37+19.37 

_3^^+a44 

-10.0 
-4.48t^2« 



77.92l}]i;°^ 
35.4ll^«.f 



_307+a62 39.74+8.60 
-2.47t»:^5t 



24.48l^:«| 



43.19 



+7.34 
-7.81 



+4.77 



102.66_^_^ 
105.94lt^o 

91.76ttll 
70.82t^:g? 
66.621^:^-^ 
45.391:111 
39.19l3:|i 
27.931^:^1 



1.03(67) 

1.10(76) 
0.98(67) 
0.93(69) 
0.92(134) 
0.88(147) 
1.16(178) 
1.28(182) 
1.24(180) 
1.09(175) 
1.23(164) 
1.26(149) 
0.95(141) 
1.14(127) 



-34- 



Table 3: Xred different models for time-resolved spectral analysis of GRB 081221 



Method 


3 Second time bins 


1 Second time bins 


{xlJ of full GRB 


ixled) (2nd pulse) 


ixlJ of full GRB 


ixled) (2nd pulse) 


BBPL (r free) 


1.31 ± 0.35 


1.52 ±0.32 


1.11 ±0.26 


1.21 ±0.28 


BBPL (r frozen) 


1.30 ±0.35 


1.50 ±0.33 


1.16 ±0.27 


1.22 ±0.29 


Band 


1.09 ±0.14 


1.17±0.11 


1.00 ±0.16 


1.04 ±0.18 


mBBPL 


1.15 ±0.14 


1.23 ±0.13 


1.07 ±0.17 


1.06 ±0.19 


2BBPL 


1.09 ±0.15 


1.17 ± 0.13 


1.02 ±0.16 


1.05 ± 0.17 



Table 4: Study of spectral evolution in pulse 2 (17.0 to 40.55 s) of GRB Q81221 (neglecting 30-40 keV) 



Model 


(dof) 




u a P 


-^peak 


P 


r 


kTh/mjkT » 


kTi » 


Case I: Count per time bin > 3000 — 


Rising part (17.0 to 21.45 s 3 bins) 












Band 


364.41 (354) 


1.0 ±0.3 


-0.441°:°^ -7.611^1^ 


98.5911^5 


— 


— 


— 


— 


BBPL 


455.92 (354) 


0.5 ±0.3 


3.2 ±1.0 — — 


— 


— 


i.89i°:gl 


24 1fi+0-64 


— 


mBBPL 


355.68 (353) 


0.9 ±0.2 


0.8 ±1.0 — — 


— 


r. o-i +0.08 
U-O-'--0.04 




40.02t|«I 


— 


2BBPL 


351.78 (352) 


0.6 ±0.1 


1.7 ± 1.0 






1 n 1+0.16 


28 7'l+i'P'7 


q 7r,+l.ll 


Case I: Count per time bin > 3000 — 


Falling part (21.55 to 40.55 s 9 bins) 












Band 


1188.17 (993) 


-2.1 ±0.1 


— -0.68 ± 0.05 -3.551^:4^ 


115.5;i;° 


— 


— 


— 


— 


BBPL 


1557.25 (993) 


-1.9 ±0.1 


-3.2 ±0.4 — — 


— 


— 


9 n9+0-03 
^•^^-0.02 


26.8lt°;t? 


— 


mBBPL 


1223.39 (992) 


-2.0 ±0.2 


3.5 ±0.3 — — 


— 


n 74+0-02 

^•'^-0.03 


9 nq+0.09 

Z.UO_Q Qg 


49.931?:!^ 


— 


2BBPL 


1147.53 (991) 


-1.9 ±0.1 


-3.1 ±0.4 — — 


— 


— 


9 1 cr+0.09 
^•^^-0.08 


38.13li:i 


iq qq+0.71 

10.00_Q gg 


Case IL Count per time bin > 1000 — 


- Rising part (17.0 to 21.45 s 10 bins) 












Band 


1247.91 (1187) 


1.5 ±0.3 


-0.44in^ -9.151^^°^ 




— 





— 


— 


BBPL 


1328.50 (1187) 


1.0 ±0.3 


4.4 ±0.8 — — 


— 


— 


1 8Q+0-05 


99 c;c+0.64 
ZZ.OO_Q g3 


— 


mBBPL 


1239.20 (1186) 


0.9 ±0.1 


1.3 ±0.5 — — 


— 


0.87 ± 0.09 


1 q6+°-2^ 

^•^"-0.14 


39.0lt?:^« 


— 


2BBPL 


1222.76 (1185) 


0.4 ±0.6 


9.0 ±2.0 — — 


— 


— 


9 Q7+0.37 
-0.16 




q ofi+l Ol 

».»U_Q gg 


Case IL 


Count per time bin > 1000 — 


- Falling part (21.55 to 40.55 s 29 bins) 












Band 


3743.36 (3448) 


-2.5 ±0.1 


— U./O_Q05 — O.OD_o.77 


125.7 ±3.9 


— 


— 


— 





BBPL 


4133.63 (3448) 


-2.1 ±0.1 


-3.6 ±0.2 — — 






2.04 ± 0.03 


28.15t°:^« 




mBBPL 


3804.05 (3447) 


-2.0 ±0.1 


3.0 ±0.5 — — 




n 79+0.03 
"■'-^-0.02 


9 99+0.11 
"^■"^^-0.08 


53.0011^^ 




2BBPL 


3690.21 (3446) 


-2.0 ±0.2 


-3.6 ±0.4 — — 






9 qi+0.13 

^■''-'--0.11 




1 "i 79+0.65 



" The values quoted are for the first time bin. 



Table 5: Study of spectral evolution in pulse 1 (-1.0 to 12.05 s) of GRB 081221. The bins are obtained by requiring > 

1000 counts per bin (30-40 keV bins are neglected) 



Model 


(dof) 


M 


u 


a 


/3 


TP °' 


p 


r 


kTh/kTin/kT ^ 


kTi » 


Rising part (-1 to 2.15 s 


1 bin) 


















Band 


115.78 (116) 






-0 

U.OO_Q 22 


-10.0 


170.31^°:? 










BBPL 


109.67 (116) 














1 gQ+O-35 


oo O7+4.08 
D0.Z/_3 7g 




mBBPL 


110.27 (115) 












98+°° 

U.y5_Q 28 




62.78l^«2^^ 




2BBPL 


103.05 (114) 














-^•'^-3.04 


38.47lt« 


6.57l?:^I 


Falling part (2.25 to 12.05 s 4 bins) 


Band 


544.65 (473) 


-0.7 ±0.1 




-0.861°:?^ 


-3.6lir 


82.021^:^1 










BBPL 


571.47 (473) 


-0.7 ±0.1 


-0.2 ±0.4 










2.09l°:l? 


19.62lJ:^^ 




mBBPL 


548.22 (472) 


-0.7 ±0.2 


2.5 ±0.8 










1.5lir 


41.6311^,^ 




2BBPL 


544.79 (471) 


-0.7 ±0.2 


-0.1 ±0.5 










2.0410:?! 


28.5911:99 


9 95+2-30 



" The values quoted are for the first time bin. 
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Table 6: Comparison of different model fits at different episodes of GRB 08 f 221 



Region 


iViOQei2 / ivioueii 


P 


(J 


P T 


Pulse 1, Rising part 


BBPL/Band 


0.385 


0.87 


61.46% 


mBBPL/Band 


0.415 


0.81 


58.50% 


2BBPL/Band 


0.298 


1.04 


70.17% 


2BBPL/BBPL 


0.029 


2.19 


97.12% 


Pulse 1, Falling part 


Band/BBPL 


0.301 


1.03 


69.93% 


mBBPL/BBPL 


0.334 


0.965 


66.57% 


2BBPL/BBPL 


1.29x10 ^ 


4.36 


99.99% 


Band/2BBPL 


0.480 


0.705 


51.95% 


Pulse 2, Rising part 
(> 3000 counts/bin) 


Band/BBPL 


0.018 


2.37 


98.23% 


mBBPL/BBPL 


0.011 


2.55 


98.93% 


2BBPL/BBPL 


1.5x10 


9.29 


100% 


2BBPL/Band 


0.390 


0.86 


60.94% 


Pulse 2, Falling part 
(> 3000 counts/bin) 


Band/BBPL 


1.04x10"^ 


4.41 


99.99% 


mBBPL/BBPL 


7.86x10"^ 


3.95 


99.99% 


2BBPL BBPL 


1.99x10 


17.22 




2BBPL/Band 


0.303 


1.03 


69.71% 


Pulse 2, Rising part 
(> 1000 counts/bin) 


Band/BBPL 


0.172 


1.36 


82.76% 


mBBPL/BBPL 


0.118 


1.56 


88.17% 


2BBPL/BBPL 


4.55x10 


9.66 


100% 


2BBPL/Band 


0.374 


0.89 


62.60% 


Pulse 2, Falling part 
(> 1000 counts/bin) 


Band/BBPL 


0.0018 


3.12 


99.82% 


mBBPL/BBPL 


0.0075 


2.67 


99.24% 


2BBPL/BBPL 


1.23x10-^5 


19.61 


100% 


2BBPL/Band 


0.343 


0.95 


65.64% 



Table 7: Study of spectral evolution in pulse 1 (-1.0 to 40.85 s) of GRB 090618. The bins are obtained by requiring > 

1000 counts per bin 



Model 


(dof) 


M 


u 


a 


/3 


a 


P 


r 


kn/kTijkT " 


kTi » 


Rising part (-1 to 14.15 


s 10 bins) 


















Band 


623.94 (547) 


-0.8 ±0.2 




-0 46+°-^° 


-3.07l°:i? 


344.9ti:J 










BBPL 


661.77 (547) 


-0.6 ±0.2 


-1.4 ±0.4 










-■ 7-1+0.05 
-■-•'^--0.04 


64.79111^ 




mBBPL 


621.93 (546) 


-0.7 ±0.2 


-0.5 ±2.0 








0.83lH^ 


■'-•'-'^-0.10 






2BBPL 


624.55 (545) 


-0.7 ±0.2 


-1.4 ±0.8 










-1 79+0.11 

-■-•'^-o.os 




25.8811^^ 


Falling part (14.15 to 40.85 s 11 bins) 


Band 


602.53 (571) 


-1.0 ±0.3 




'-'•'^-0.11 




186.2ti«| 










BBPL 


642.37 (571) 


-0.9 ±0.3 


-1.4 ±0.7 










1.83 ± 0.05 


37.95;2J7 




mBBPL 


602.30 (570) 


-1.1 ±0.2 


1.5 ± 1.5 








69+° °^ 


-■■•'^-0.19 


85.8ll^^gf 




2BBPL 


593.74 (569) 


-0.3 ±0.3 


-2.9 ± 1.0 










o 1 1+0.10 
^•-•-^-0.12 


50.88lt«8 


1 ^ 00+1.93 



" The values quoted are for the first time bin. 



Table 8: Study of spectral evolution in 


pulse 2 (61 to 75.0 s) of GRB 090618. 


The bins 


are obtained by requiring > 




2000 counts per 


bin 














Model (dof) 




u 




P 


r 


kn/kTin/kT " 


kTi « 


Rising part (61 to 64.35 s 1 bin) 


Band 2012.93 (1707) 


8.0 ±3.0 




-0.68lH^ -2.49in^ 192.5ilJJ 










BBPL 2173.47 (1707) 


5.0 ±1.0 


6.5 ±5.0 






1.64 ±0.02 


42.82^J:}^ 




niBBPL 2099.21 (1706) 


8.0 ± 1.5 


10.0 ±5.0 




0.68 ± 0.02 


1.29 ± 0.08 


97.1;!:^ 




2BBPL 2030.83 (1705) 


6.0 ± 1.5 


3.0 ±3.0 






-1 70+0.05 

-'-'°-0.04 


lie; c;(?+13.33 


33.19li;^i 


Falling part (64.35 to 74.95 


s 24 bins) 














Band 1574.79 (1317) 


-14.0 ± 1.0 




-0.88 ± 0.03 -2.74toil 262.8^^1 










BBPL 2402.54 (1317) 


-5.3 ±0.3 


-6.9 ±0.9 






1.78 ± 0.01 


p-o 01 +1.01 

oo-y-L-i.oo 




mBBPL 1794.05 (1316) 


-8.0 ±3.0 


15.0 ±5.0 




•"••"^-O.OOS 


1 f^8+0-03 
^■^°-o.io 


157.36lt^^ 




2BBPL 1794.06 (1315) 


-5.5 ±0.5 


-4.0 ±0.5 






1.78 ± 0.02 




01 70+1.68 



" The values quoted are for the first time bin. 
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Fig. 1. — Background subtracted Light curve (LC) of GRB 081221, fitted with Norris model 
(Norris et al. 2005). The LC is generated by adding the two highest count Nal (nl and n2) 
and one EGO (bO) detectors after binning by 2 s. 
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Fig. 2. — Comparison of Xred of Band, BBPL, mBBPL and 2BBPL in the time-resolved 
spectroscopy of GRB 081221. The symbols are explained in the inset. For convenience, we 
draw line to join the BBPL (F free) and dot-dashed line to join 2BBPL points. 
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Fig. 3. — Time evolution of kT and Epeak of GRB 081221. The filled circles are the BBPL 
with r free, while the open circles are BBPL with F frozen. The pluses are Epeak- The 
light curve (LC) of the GRB is shown as a histogram with errors. We have divided the LC 
into three regions. We see a clear hard-to-soft evolution in pulse 1. The same evolution 
is seen in the falling part of the second pulse. In the overlapping region, the variation is 
rather soft-to-hard. This might be the effect of overlapping, or the evolution may as well be 
genuinely intensity tracking (see text). 
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Fig. 4.— Verification of Liang & Kargatis (1996, LK96) law in tlie pulses of GRB 081221. 
The first pulse (left panels) shows a strict hard-to-soft evolution, while the second pulse is 
intensity tracking. The x-axis represents the "running fluence", defined by LK96. For the 
BBPL model, fluence here means that of the BE component. 
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Fig. 5. — Flux evolution of GRB 081221. Crosses represent BB flux, while triangles represent 
PL flux. The total flux is marked by open boxes. Energy flux, in the units of 10~^ erg cm~^ 
s~^, is plotted in the upper panels; photon flux, in the units of photon cm~^ s"^, is plotted 
in the lower panel, for both T free and frozen cases. 
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Fig. 6. — Evolution of a (triangles), (3 (stars) of Band and F (pluses) of BBPL throughout 
the GRB. The error bars of parameter (5, being large and undetermined in many cases, are 
not shown. In some cases, the /3 pegs at -10.0. Notice that the parameters are more or less 
constant throughout the burst. Hence, they can be tied (see text). 
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Fig. 7.— Fitting the 29.0 to 32.0 s spectral data of GRB 081221 with BBPL (lower panels) 
and 2BBPL (upper panels) models with the 30-40 keV band included (left panels) and 
neglected (right panels). uF^ has the unit keV^( Photon cm~^ s~^). Note the double hump 
structure in the residuals of the BBPL model. This structure is taken care of by the 2BBPL 
model. 
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Fig. 8. — Predicted peak energy {Ep^ak) of Band, kT of BBPL, kTin of mBBPL and kT^ 
of 2BBPL models are compared with the observed values. The open circles are the values 
obtained for wider bin size (> 3000 counts per bin), while the filled circles represent values 
obtained for finer bin size (> 1000 counts per bin). The dot-dashed hne is the line of equality. 



-48 - 




100 



fc^^ 10 




2 

X 



-2 - 



10 



100 



1000 



1 


h-^ — \ — \ — h-h 


1 1 1 1 1 1 

^^^^ 




1 1 III 








-HHH \ H- 




MM 


. 1 


H 


. . . 1 




,1 , , , 





10 



100 



1000 



Energy (keV) 



Fig. 9. — Comparison of spectral fitting between 2BBPL (upper panels) and Band (lower 
panels) model for -1.0 to 2.15 s time bin of GRB 081221. uF^ has the unit keV^ (Photon 
cm~^ s""^). The right panels show the fit with the 30-40 keV channels neglected. Note the 
structure in the residual of the Band model — positive excess near 15 keV and 150 keV, 
and negative excess near 40-60 keV. Compared to this the 2BBPL model does not show 
any structures in the residual. The 2BBPL model is preferred over the Band model at 1.04 
a with 70% confidence level and p-value = 0.298 for 30-40 keV neglected case and at 0.95 
a with 65.5% confidence level and p-value = 0.341 for 30-40 keV included case, based on 
F-tests. 
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Fig. 10. — The residual in various time bins of the falhng part of GRB 081221. The data 
used are those with > 3000 counts per bin. The time bins used are 2nd, 6th and 9th bins 
(from top to bottom). The residuals are obtained by omiting the lower BB from the 2BBPL 
fit. The lower BB models are overplotted with the residuals to show the significance of this 
BB component. Different lines and symbols show different detectors. Note that the lower 
BB temperature shifts to the lower energy with time. 



